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Abstract 

Quantum fluctuation is introduced into the Markov random fields (MRF's) 
model for image restoration in the context of Bayesian approach. We in- 
vestigate the dependence of the quantum fluctuation on the quality of BW 
image restoration by making use of statistical mechanics. We find that the 
maximum posterior marginal (MPM) estimate based on the quantum fluc- 
tuation gives a fine restoration in comparison with the maximum a posterior 
(MAP) estimate or the thermal fluctuation based MPM estimate. 
PACS numbers : 02.50.-r, 05.20.-y, 05.50.-q 
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I. INTRODUCTION 



Recently, the problems of information science were investigated from statistical me- 
chanical point of view. Among them, the image restoration is one of the most suitable 
subjects. In the standard approach to the image restoration, an estimate of the original 
image is given by maximizing a posterior probability distribution (the MAP estimate) |]T[| . 
In the context of statistical mechanics, this approach corresponds to finding the ground 
state configuration of the effective Hamiltonian for some spin system under the random 
fields. On the other hand, it is possible to construct another strategy to infer the orig- 
inal image using the thermal equilibrium state of the Hamiltonian. From the Bayesian 
statistical point of view, the finite temperature restoration coincides with maximizing a 
posterior marginal distribution (the MPM estimate 0,0) and using this strategy, the er- 
ror for each pixel may become smaller than that of the MAP estimate. As we use the 
average of each pixel (spin) over the Boltzmann-Gibbs distribution at a specific tempera- 
ture, the thermal fluctuation should play an important role in the MPM estimate. Then, 
the temperature controls the shape of the distribution and if we choose the temperature 
appropriately, the sampling from the distribution generates the important configurations 
for a fine restoration. Besides this hill-climbing mechanism by the thermal fluctuation, 
we may use another type of fluctuation, namely, the quantum fluctuation which leads to 
quantum tunneling between the states. If we use the sampling from the Boltzmann-Gibbs 
distribution based on the quantum fluctuation, it may be possible to obtain much more 
effective configurations for a good restoration. The idea of the MRF's model using the 
quantum fluctuation was recently proposed by Tanaka and Horiguchi [Q, however, they 
investigated the quantum fluctuation in the context of the optimization (the MAP es- 
timate by the quantum fluctuation) and they used the ground state as the estimate of 
the original image. We would like to stress that we use the distribution based on the 
quantum fluctuation itself and the expectation value is used to infer the original image. 
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It is highly non-trivial problem to investigate whether the MPM estimate based on the 
quantum fluctuation becomes better than the MAP estimate or the thermal fluctuation 
based MPM estimate. 

This is a basic concept of this paper. This paper is organized as follows. In the next 
Sec. II, we explain our model system and the basic idea of our method in detail. In 
Sec. II, we also introduce the criterion of the restoration, that is, the overlap between 
the original image and the result of the restoration. In Sec. Ill, we introduce the inflnite 
range model in order to obtain analytical results on the performance of the restoration, 
and calculate the overlap explicitly. In Sec. IV, we show that quantum Monte Carlo 
simulations in 2-dimension support our analytical results. In Sec. V, we introduce the 
iterative algorithm which is derived by mean-fleld approximation and apply this algorithm 
to image restoration for standard pictures. The last Sec. VI is devoted to discussion about 
all results we obtain. In this section, we also mention the inequality which gives the upper 
bound of the overlap. 



II. BASIC IDEA AND FORMULATION 

Let us suppose that the original image is represented by a conflguration of Ising spins 
{^} = = ±l;i = l,---,iV} with probability Ps({0). These images are sent 

through the noisy channel by the form of sequence {^}. Then, we regard the output of 
the sequence {^} through the noisy channel as {r}. The output probability for the case 
of the binary symmetric channel (BSC) is specified by the following form; 

We easily understand the relevance of this expression for the BSC; Lets suppose that 
each pixel changes its sign with probability and remains with 1 — Pr during the 
transmission, that is. 
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P(.. = -«.ia=P,-^^ (2) 

P(r. = 6lf.) = l-P.-^. (3) 

We easily see that there is a simple relation between flip probability Pr and inverse 
temperature as exp(2/3^) = (1 — Pr)/pT- This is reason why we refer to this type 
of noise as binary symmetric channel. Using the assumption that each pixel C,i in the 
original image {^} is corrupted independently (so-called memory-less channel), namely, 
= Hi PiTMi): W6 obtain Eq. ([1|). This BSC is simply extended to the follow- 
ing Gaussian channel (GC) 

nut({r}|{0) = ^^^^^ exp (-^E(^^ -^0^^)') • (4) 

where r is a standard deviation of observable (corrupted pixel) from scaled original 
pixel To^i. 

Then, the posterior probability P({o"}|{r}), which is the probability that the source 
sequence is {a} provided that the output is {r}, leads to 

P(^.^\^^^^ Pi{r}\{a})PmiW}) ... 

J2aP{{r}\{(T})PUW}) 
by the Bayes theorem. As we treat the BW image and the BSC (0), a likelihood 
P({r}|{,^}) is appropriately written by 



P({r}|{0)~exp (^hJ2r.a}j 



(6) 



PmiW}) appearing in the Bayesian formula (^) is a model of the prior distribution Ps({^}) 
and we usually use the following type; 

PmiW}) ~ exp (3m ^i^j C^) 

\ <ij> J 

where J2<ij>{- " ") means the sum with respect to the nearest neighboring pixels and 
controls the smoothness of the picture according to our assumption. Substituting Eqs. 
(I) and (p2D into Eq. (P), we obtain the posterior probability P({o"}|{t}) explicitly; 
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exp E<ij> o-iaj + hJ2i no-i 
In the framework of the MAP estimate, we regard a configuration {a} which maximizes 
the posterior probabihty -P({o"}|{t}) as an estimate of the original image {^}. Obviously, 
this estimate {a} corresponds to the ground state of the following effective Hamiltonian 
(the random field Ising model) 

<ij> i 

Therefore, in the limit of f3rn/ h — > oo, we expect that the original image should be complete 
BLACK picture or complete WHITE picture, whereas in the limit of (3^/ h^O,we assume 
that the original image should be identical to the observable {r} itself. 

On the other hand, in the framework of the MPM estimate, we maximize the following 
posterior marginal probability 

p/ ij ^^ E.^.,P({r}|{a})P^({a}) 
EaP{{^}\W})Pm{{(r}) 

As we treat the case of BW image, the estimate of the i-th pixel should be given as 



^ ? E.^.P(MIW)Pn.({a}) \ 

^ I j:.p{{r}\w})RUW}) J 

= ^M{(^i)h,f3j (11) 

where (■ ■ ■)h,f3m means the average over the posterior probability Eq. (H). Consequently, 
our problem is reduced to that of statistical mechanics which is described by the effective 
Hamiltonian Eq. (|^). As the Hamiltonian Eq. (^) has lots of local minima due to the 
quenched disorder {r}, in general, it is quite difficult to obtain the thermal equilibrium 
state which contributes to fine restoration without being trapped in a local minimum for 
a long time. In order to overcome this difficulty, we add the quantum transverse field 

-r^af^Tii (12) 



to the effective Hamiltonian (|^) as quantum fluctuation. In this expression, af means the 
x-component of the Pauh matrix and F controls the width of the quantum fluctuation. 
Intuitively, the term F af is regarded as the tunneling probability between the eigenstates 
of the operator (2;-component of the Pauli matrix), namely, \cr^ = ±1 >. The tunneling 
probability between the states |(jf = ±1 > leads to | < erf = +1 1 T af \ af = — 1 > p = 
F^. As the result, the term (|l^) generates the superposition of the states \af = +1 > 
(BLACK) and \af = — 1 > (WHITE). Using this fuzzy representation for each pixel, 
we may construct the algorithm which is robust for the choice of the hyper-parameters, 
especially, for the edge parts of a given picture. 

Our problem is now reduced to that of quantum statistical mechanics for the next 
effective Hamiltonian 

Tieff = E -PmY. - T af = + Hi (13) 

i <ij> i 

where we defined Hi = Hcs — Hq. Our main goal is to calculate the local magnetization 
{d'f)h,i3m,r of the system described by the above Hamiltonian, that is to say, 

{(^i }h,f3„.r = — ^ (14) 

Tr<^exp(-7YefT) 

and regard the quantity sgn{{af)h,f3^,r) as an estimate of the original pixel ^j. Therefore, 
the averaged performance of our method is measured by the following overlap M{h, f3mi F) 
as 

M(/.,/?^,F) = Tr|5,.|P,({e})Pout({r}|{e})e.sgn((af),,^„,r). (15) 
Then, our main interests are summarized as follows. 

• Is it possible for us to use the quantum fluctuation in place of the thermal one ? 

• Does there exist a specific choice of F which gives the optimal image restoration ? 
Before we calculate the above overlap ([T5|), we may add the parity check term, which was 



recently introduced by Nishimori and Wong |p, to the effective Hamiltonian (|^). This 
parity check term is represented as PjJ2<ij> Jij^t^j^ a^id we rewrite Tio as 
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^0 = E J^j^t^] - E ^.''^I - ^ E ^^^t (16) 

<ij> <ii> i 

where Jij is the noisy version of the product of arbitrary two original pixels ^j^j and the 
output of this quantity through the noisy channel is given by 

^out({J}|{0) = (2eosh/?.)^^ (^^' ^ ^--^^^^-^ (^^^ 

for the BSC and 



<ij> 



Pont{{J}m) = ^^j^^, exp (^-^ EJJ,, - Joe.e,^^ 

for the GC, respectively. Nb is the number of the terms appearing in the sum in Eq. 
(|T7|) or Eq. (|l^). Then, the effective Hamiltonian Hcs = Ti-o + Ti-i describes the thermo- 
dynamics of quantum spin glass under random fields. 

In the next section, we introduce the rather artificial model, namely, the infinite range 
model in which spins in the system (0) are fully connected. 



III. THE INFINITE RANGE MODEL 

In this section, we calculate the overlap (|T^) explicitly using the infinite range version 
of the effective Hamiltonian (|T^). We use the GC for the analysis of the infinite range 
model in this section and the the BSC for the quantum Monte Carlo simulations in the 
next Sec. IV, respectively. However, these two channels can be treated by the following 
single form. 

^out({J}|{r}) = n FMj) n ^l(^..)exp [f3r E J^MJ+PrT.'^^^i] (19) 
<ij> <ij> \ <ij> i I 



with 
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for the BSC and with 
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i^i(n) = ^exp(-^(rf + ro^)), (21) 

for the GC, and we set (3j = Jo/J'^,Pr = tq/t^- 

As the original image, we use the ferro-magnetic snapshot from the distribution 

P.(K})^^exp (22) 
where ■ ■) means the sum over all possible combinations of and we divided the 



argument of the exponential in Eq. (22) by to take a proper thermo-dynamical limit 
as Hamiltonian should be of order A^. For the same reason, we should re-scale the terms 
appearing in Eq. (|13|) as PjY.<ij> JijC^t^j iPj/N)Y.ijJij^!cr^ and PniJ2<ij> cr^cr^j 
{13m /N) cr-aj when we treat the infinite range model. 

It must be noted that Tio and Ti.i do not commute with each other and we use the 
following the Trotter decomposition 

Z= limTr^. e" — e"— (23) 

P— »oo \ / 

to calculate the partition function explicitly. In this formula, Tio and Tii are eigenvalues 
of the operators TYq and Tii with respect to the following eigenvector 

N 

\Wl}>=X{®\<> ik = l,.-;P) (24) 

i=l 

with 

^tk\(^!k > = ^ik^ik > ■ (25) 

P means the Trotter number and we distinguish the different Trotter slices by the indices 
k. 

Now we can calculate the partition function for the quantum spin system (^) in terms 
of the corresponding classical spin system whose dimension increases by 1. Using the 
Trotter formula (the path integral formula) and well-known replica method [Q, namely. 



lim 

n— >oo 



1 



n 



(26) 



we can obtain the overlap as a function of the macroscopic parameters j3m and F by 
making use of the saddle point method. The bracket [■ ■ ■] denotes the average over the 
distribution P,({0)Pout({ J}, {r}m)- 

The standard replica calculations and saddle point method lead to the following cou- 
pled equations. 



[ii] =mo= tanh(/3omo) 



[6('^iii')fe,/3m,r] - t 



^-^TT-o -J Dun~' Duj^y-'smhy 

2 cosh(psmo) J-oo J-oo 



-oo 



' — oo 
coo 



2 cosh(/5smo) J-oo 

Tr^ e^'"^"^ 
2 cosh.{j3smo) J- 



/oo roo 
-oo J —oo 



oo 
oo 



' —oo 
roo 



/oo r poo 
Du il'^ / Duj^y^^sinhy 
-oo L J—oo 



2 cosh(/5smo) J-oo 



-oo 
roo 



(27) 
(28) 

(29) 

(30) 



/oo /-oo 
Dwr^-i / Du <^^y-\oshy 
-oo J—oo 



+ r 



/oo 
Z^cj y^^sinhy 
-oo 



(31) 



where (■ ■ ■)h,f3m,T means the average by the posterior probability using the same way as 
Eq. (|n|). Du or Dy means Gaussian integral measure Du = due"^^^"^ /\/2n. In order 
to obtain the above saddle point equations, we used the replica symmetric and the static 
approximation, that is. 



tK = t 

Sa{KL) = S{K ^ L), 1{K = L) 

QaP = Q- 



We also defined functions $, y and Q as 



$ = u^irhy + Q{J(3.jY + J(3uj^S - g + (ro/i + J^f3M + 



/oo 
Duj coshy. 
-oo 



(32) 
(33) 
(34) 

(35) 
(36) 
(37) 



Then the overlap which is a measure of retrieval quality is calculated explicitly as 

Trc f 6^""*°^ r°° r°° 



xsgn 



u^{Th)^ + Q{jPjY + (ro/i + JoPM + P^rn + J^jw^S-Q 



(St 



where the above overlap M depends on F through m (pSl). 

We first consider the case of (3j = 0, that is to say, the conventional image restora- 
tion. We choose a snapshot from the distribution (p2|) at source temperature = 0.9. 
According to Nishimori and Wong [H, we fix the ratio h/Pm and adjust /3m(= as 
a parameter for simulated annealing |10| and controls F as a quantum fluctuation. If we 



set F = 0, the lines of M{Tm,T = 0) should be identical to the results by the thermal 
MPM estimate 0. On the other hand, if we choose = and F = 0, the resultant 
line M{Tm = 0, F) represents the performance of the quantum MAP estimate. We should 
draw attention to the fact that the quantum fluctuation vanishes at F = 0. In practical 
applications of the quantum annealing based on quantum Monte Carlo simulations, 
we should reduce F from F > to F = during Monte Carlo updates. However, the 
resultant performance obtained here is calculated analytically provided that the system 
reaches its equilibrium state. Therefore, we can regard the result M{Tm = 0, F = 0) as a 
performance when F is decreased slowly enough. 

In FIG. ^ we set the ratio h/ [3m to its optimal value Pr/Ps = 0.9 and plot the overlap 
M{Tm, F) for the case of F = 0, 0.5 and F = 1.0. Obviously, for the case of F = 0, the 
maximum is obtained at a specific temperature = 0.9(= Tg) 0. However, if we add a 
finite quantum fluctuation, the optimal temperature is shifted to the low temperature 
region. 

In FIG. 0, we plot M{Tm, F) for the case of = 0, 0.1, 0.9 with the fixed optimal ratio 
h/Pm = 0.9. This figure shows that if we set the parameters h, Pm to their optimal value 
in the thermal MPM estimate, the quantum fluctuation added to the system destroys 
the recovered image (see the lines M(Tm,,F) for the case of Tm = 0.9). Therefore, we 



may say that it is impossible to choose all parameters h, jSm and F so as to obtain the 
overlap which is larger than Mmax = MiT^ = 0.9, F = 0). This fact is also shown by 
3-dimensional plot M{Tm, F) in FIG. 0. 

Although, we found that a finite F does not give the absolute maximum of the overlap, 
the quantum MPM estimate M{Tm = 0, F > 0) has another kind of advantages. As FIG. 
^ indicates, the overlap of the the quantum MPM estimate is almost fiat in comparison 
with M(Tm = 0.1, F > 0) or M{Tm = 0.9, F > 0). This is a desirable property from 
practical point of view. This is because the estimation of the hyper-parameters is one of 
the crucial problems to infer the original image, and in general, it is difficult to estimate 
them beforehand. Therefore, this robustness for hyper-parameter selection is a desirable 
property. We also see this property in FIG. |^. 

As we already mentioned, the overlap at = and F = corresponds to the result 
which is obtained by quantum annealing |]12[, that is to say, the quantum MAP estimate. 
We see that the result of the quantum MPM estimate is slightly better than that of the 
quantum MAP estimate. 

We next show the effect of the parity check term. In FIG. we set = Tg = 0.9, h = 
1.0 and Jq = J = 1.0 and plot the overlap as a function of /3j for several values of F. We 
see that the performance of the restoration is improved by introducing the parity check 
term which has much information about the local structure of the original image. 

In the next section, we check the usefulness of this method in terms of quantum Monte 
Carlo simulation. 



IV. QUANTUM MONTE CARLO SIMULATION 

In this section, Monte Carlo simulations in realistic 2-dimension are carried out in 
order to check the practical usefulness of our method. We use the standard pictures which 



are provided on the web site [[I^] as the original image, instead of the Ising snapshots. 



In order to sampling the important points which contribute to the local magnetization 
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(af), we use the quantum Monte Carlo method which was proposed by Suzuki [Q. As we 
mentioned in the previous sections, we can treat the (i-dimensional quantum system as 
{d + l)-dimensional classical system by the Trotter decomposition In this sense, the 
transition probability of the Metropolis algorithm leads to 



P(a — s> <T ) = min 



l,exp{-iEi(T')-Eia))) (39) 



where E{(t) is energy of the classical spin system in {d + l)-dimension (in the present 
case, (2 + 1) = 3-dimension) as follows. 

ijk 

h 

— -p X! '^i,j^i,3,k ~ B^Yl '^i,j,k(^id,k+i (40) 
ijk ijk 



where we defined B = lncosh(r/P). The transition probability Eq. (PP|) with Eq. 



generates the Boltzmann-Gibbs distribution asymptotically and using the importance 
sampling from the distribution, we can calculate the expectation value of the i-th spin 
af , namely, {d'l)h,(3^,v^ and using this result, we obtain an estimate of the i-th. pixel of 
the original image as sgn(((3"f )/i,/3^j). We show the results in FIG.s ^ and |[ From these 
Figures, we see that there exists the optimal value of the transverse field F . In FIG.s 
1^ and ^, we display the results by quantum Monte Carlo simulations when we add the 
parity check term for the parameter sets F = 2.0, /i = 1.0 and (3m = 0.5. We see that the 
resultant pictures using the parity check term are almost perfect (see (3j = 1.0 and 1.5). 



V. MEAN-FIELD ALGORITHM 

In the previous sections, we see that the quantum fluctuation works effectively on 
image restoration problems in the sense that the quantum fluctuation suppress the error 
of the hyper-parameter's estimation in the Markov random fields model. In addition, 
by making use of the quantum Monte Carlo simulations, we could apply it to the image 
restoration of the 2-dimensional standard pictures. However, in order to carry out the 
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simulations, it takes quite long time to obtain the average {o'f)h,i3m,T and it is not suitable 
for practical situations. 

In this section, in order to overcome this computational time intractability, we derive 
the iterative algorithm based on the mean-field approximation. This algorithm shows fast 
convergence to the approximate solution. 

Within the mean- field approximation, we rewrite the density matrix p = e~^''^ /Z for 
2-dimensional version of the effective Hamiltonian TYeff as 



p ~ l[<S>pij: 



(41) 



where we defined pij as 



n=l 



(42) 



with 



-\ij{2) 



(43) 



In the above expressions, \ij{n), n = 1,2 means eigenvalues of the 2x2 matrix Hij 
[Hij]u = H^^\ [Hij]22 = Hii\ [Hij]i2 = [Hij]2i = -r) and H^f is defined by 



(t) 



(44) 
(45) 



Using this decoupled density matrix, the local magnetization at a site (i, j), namely, m^^ 
leads to 



+ 



2cosh(VQ;2 + r2) 
2cosh(7a^^r^ 



q; 



+ + P)2 - 



- Vq;2 + r2)2 + p 



(46) 



12 



For this local magnetization (^Bf), the estimate of the pixel is obtained as sgn[mjj]. 
We solve the mean-field equations (|46|) with respect to niij until the condition 

e,j = \m^^^ -mfj>\ < 10"^ (47) 

holds for all pixels {i,j}- We show its performance in FIG. ^ and TABLE |. From 
TABLE I, we see that if we introduce appropriate quantum fiuctuation, the performance 
is remarkably improved, and in addition, the speed of the convergence becomes much 
faster. However, if we add the quantum fiuctuation too much, the fiuctuation destroys 
the recovered image. We also see that the optimal value of F exists around F ~ 1.6. 



VI. SUMMARY AND DISCUSSION 

In this paper, we investigated to what extent the quantum fiuctuation works effectively 
on image restoration. For this purpose, we introduced an analytically solvable model, 
that is, the infinite range version of the MRF's model. We applied the technique of 
statistical mechanics to this model and derived the overlap explicitly. We found that the 
quantum fiuctuation improves the quality of the image restoration dramatically at a low 
temperature region. In this sense, the error of the estimation for the hyper-parameters 
Pm, h can be suppressed by the quantum fiuctuation. 

However, we also found that the maximum value of the overlap M'^^^^^^™'> never 
exceeds that of the classical Ising case may show this fact by the following 

arguments; First of all, the upper bound of the overlap for the classical system is given 
by setting h = P^-, Pg = P^, that is, 

Mit''"*''(/3r,ft) = Tl'M6e'''S.--«-P(tt})sgn|Tr„CTie'''E. p„({„})] 

|Tr<,aie^- 2..^«'^«P^({a})| 
= TTr\TT^a,e^^^^^^"^Pm{W})\. (48) 

For the quantum system, the overlap is bounded by this maximum value M^}^^^^^^^ as 
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^(quantum) _ "'«'P({^}) sgn[Tr^(j,^ e'^S, ^?P„(a^)] 

= Tr,|Tr5e.e^-S.-»«^P({0)| = Mt^r""''). (49) 
We can see this inequality more directly as follows. 

= Tr|.,^}e,e^^S^-'«»P({0)sgn[Tr^a,^e'^S.-«^?+ri:^^^p^(|^.|)] 

= M(^"'°)(/i,P^,r), (50) 

where the identity sgn(x) = x/\x\ was used. We should notice that in the left hand side 
of the above inequality (|50|), the arguments of the trace w. r. t. r always take positive 
values, while in the right hand side, they can be negative. 

In order to check the usefulness of the method, we carried out quantum Monte Carlo 
simulations in realistic 2-dimension. We found that the results by the simulation support 
qualitative behavior of the analytical expressions for overlap. 

We introduced the iterative algorithm in terms of the mean-field approximation and 
applied it to image restoration of the standard pictures. We found that the quantum 
fluctuation suppress the error of the hyper-parameter estimation. In addition, we found 
that the speed of the convergence to the solution is accelerated by the quantum fluctuation. 

From all results obtained in this paper, we concluded that the quantum fluctuation 
turns out to enhance tolerance against uncertainties in hyper-parameter estimation. How- 
ever, if much higher quantities of restoration are required, we must estimate those parame- 
ters using some methods. One of the strategies for this purpose is selecting the parameters 
f3mi h and F which maximize a marginal likelihood. By making use of the infinite range 
model, the usefulness of this method can be evaluated. The details of the analysis will be 
reported in forth coming paper. 

Of course, the application of this strategy to the restoration of gray-scaled image 
13| , p!^ will be considered as an important future problem. 
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FIGURES 

FIG. 1. The overlap M as a function of Tm for several values of F. We set the system 
parameters as = 0.9, tq = r = 1.0 and h/[5m = 0.9 = (3r/Ps- For F = case, the optimal 
temperature coincides with the source temperature Tg = 0.9. As the quantum fluctuation 
F increases, the optimal temperature is shifted to a low temperature region. However, the 
maximum value of the overlap does not change. 

FIG. 2. The overlap M as a function of F for several values of T^- We set the system 
parameters as Tg = 0.9, tq = t = 1.0 and h/Pm = 0.9 = Pr/f^s- The overlap at Tm = and 
F = corresponds to the result by quantum annealing. The quantum MPM estimate works 
effectively at a low temperature region and the results are robust for the choice of F. 

FIG. 3. The overlap M as a function of the quantum fluctuation F and temperature 

FIG. 4. The overlap Af as a function of (3j for several values of F. We set the system 
parameters Tm = Tg = 0.9, r = tq = 1.0, J = Jq = 1.0 and h/Pm = 0.9 = Pt/Ps- For the case 
of F = 0, the optimal (3j is naturally identical to Jq/ = 1.0. As the quantum fluctuation 
F increases, the overlap M decreases because the quantum fluctuation destroys the recovered 
image. 

FIG. 5. The results by quantum Monte Carlo simulations for standard picture ( A Japanese 
Kanji stamp for the name of suzuki which is the most popular name in Japan. The size is 
50 X 50.). From the upper left to the lower right, the original picture, the damaged picture, the 
results of F = 0.2,1.0,1.7,1.9,2.1 and 2.7 are displayed. The noise rate is 10% (The overlap 
between the original picture and the damaged one is 0.9). 

FIG. 6. The overlap M as a function of the quantum fluctuation F for the standard picture 
in FIG. p. We set Pm = 0.5, h = 1.0 and P = 50. The error-bars are calculated by averaging 
over five independent runs. 
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FIG. 7. The results by quantum Monte Carlo simulations including the parity check term. 
We fixed T = 2.0, h = 1.0 and = 0.5. From the upper left to the lower right, the original 
picture, the damaged picture, the results of f3j = 0.01,0.5,1.0 and 1.5 are shown. The noise 
rate is 10% (The overlap between the original picture and the damaged one is 0.9). 

FIG. 8. /3j-dependence of the overlap M calculated by quantum Monte Carlo simulations 
for standard picture in FIG. We set f3m = 0.5, /i = 1.0 and F = 2.0. The error-bars are 
calculated by five independent runs. 

FIG. 9. The restored pictures ( Their size are 50 x 50) by quantum iterative algorithm for 
several values of F. From the upper left to the lower right, the original image, the corrupted 
image, the results of F = 0.001, 0.8, 1.2, 1.6, 2.0 and F = 3.0. The noise rate is 20% (The overlap 
between the original picture and the damaged one is 0.8). 
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TABLES 

TABLE L The overlap M calculated by the quantum iterative algorithm for several values 
of r. The restored pictures are shown in FIG. ^. The iteration times are also listed. We set 
T = l,J = 0.5. 
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